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ABSTRACT 


The general differential equations describing the behavior of the 
rotary regenerator including longitudinal heat conduction in the direction 
of fluid flow are sufficiently complicated to preclude a complete analyti- 
cal solution. A few solutions found in the Literature restricted to very 
special cases are discussed, 

A numerical finite cirference method is presented which will 
determine the effect of Longitudinal heat conduction in rotary regenerators 
for steady state conditions. In the development no assumptions are made 
which would restrict the range of parameters for which the analysis 
would be applicable. 

The conduction effect on the regenerator effectiveness was evaluated 
with the view of obtaining results most useful for the gas turbine 
regenerator problem, however, these results may also be used for other 
regenerator problems. 

A CDC 1604 digital computer was used to carry out the computations. 
The results are presented graphically and in tabular form, employing a 
suitable set of non dimensional parameters. The range of parameters 


which have been covered arc: 


0.5 = Cnin/Cmax =1.0 


10.0 ت1۰8 


1.0 NTUo 420.0 
Y 
0.25 < (hA) 1.2 
0 Dum 1072 
0.15 As". 1.3 


The author expresses hís appreciation to C. P. Howard, Associate 


Professor, for his direction and encouragement in this work. 
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KOMEECUA YU E 
a = heat transfer arca on side deatguatcd by subseripi, sq ft 
AS, = Solid area available for longitudinal heat conduction on the side 
Of Cmax, 8q ft 


As. = solid area available for longitudinal heat conduction on the side 


n 
ОЕ Си, 64 ЕК 


As = total solid area available for longitudinal heat conduction, sq ft 


03 
в 


heat capacity rate (Wc) of fluid or rotor matrix according to 

subscript, Btu/hr, deg F 

c = specific heat of fluid (constant pressure) or rotor matrix material 
depending on subscript, Btu/lb, deg F 

h = unit conductance for thermal convection heat transfer, Btu/hr, 
deg F, sq ft 

k = unit thermal conductivity, Btu/hr, sq ft, deg F/ft 

L = total length of the matrix in the direction of fluid flow, ft 

N = number of subdivisions according to subscript 

Q = heat transfer rate, Btu/hr 

T = temperature of matrix, deg F 

(ii) S= subscripted matrix temperature on the side of max, deg F 

Ta(£,8) = subacripted matrix temperature on the side of Cmin, deg Г 

t = temperature of fluid, deg F 

tx(i,j) " subscripted fluid temperature on the side of Cmax, deg F 

Un(f£,g) = gubscripted fluid temperature on the side of Coin, deg F 

М = mass flow rate of fluid (lb/hr) or mtrix (Lb, rev/hr), according 
to subscript 

y = distance of a point from the cold end, measured in units of total 


length of matrix, O<y L 


vil 





UL = time measured in units of one halt period 


ΔΙ, А2,..............; А, 


Сі, С2,. e e e o ө ө ө е е . , Cy, 

01, Ә2,.............. ‚ 7ر‎ = constants. 
Е], Е2,. ° ° @ ° o e e e @ ө е > Ey 

Bp DERE e e e e ө е е 9 ο E, 

Subscripts.- 


avg = average 

c = cold side 

h = hot side 

i = inlet 

min - minimum magnitude 
max - maximum magnitude 
n - side of Cin 

o = outlet 

r = rotor 


side of Салах 


“© = subscript an e to indicate value extrapolated to an infinite 


x 


number of elements. 


Dimensionless parameters. - 
е = exchanger heat transfer effectiveness, ratio of actual to thermo- 


dynamically Limited maximum possible heat transfer rate 


= АА = = longitudinal heat conduction effect 


viii 





Сар пах = Capacity rate retio of fluid flow streams 


capacity rate ratio of rotor matrix to minimum fluid capacity‏ = ہا /مة 


rate 


(hA)" » (hA)n/(hA),, conductance ratio 


NTU = (hA)/C, number of transfer units on side designated by subscript 


NTUo = NTUn | ——~ , over-all number of transfer units 


η" conduction parameter on side of Cmax 


eg A conduction parameters on side of Cyn 


A > total conduction parameter 


As” = As,/As,, conduction area ratic, ratio of solid areas available 


for heat conduction in the direction of fluid flow 


ік 





a Introduction 

The regenerative type of heat exchanger has been subjected to 
extensive treatment vithin the last thirty years cr so and contributions 
to the theory of regenerators have been made by many writers. The most 
questionable of the idealizations made in the development of the theory 
is that of zero matrix conductivity in the direction of fluid flow. The 
loss due to heat conduction in the matrix material is specially important 
in the case of high performance regenerators in which a small reduction 
in the efiectiveness will result in a large increase of the heat transfer 
areas to compensate for the loss. Another important factor is the modern 
trend toward smaller sizes, since the loss due to conduction of heat in 
the matrix obviously increases with decreaging Length. 

Schultz P solved the problem including conduction for the 
special case of the balanced regenerator and very large values of Cr/Cmin 
(i.e., the capacity rate ratio of rotor matrix to minimum fluid capacity 
rate). In this case the system of partial differential cquations defining 
the problem for the steady state are reduced to ordinary differential 
equations and the problem can be solved analytically. 

Hahnemann ЕЗ extended the solution to the unbalanced regenerator 
(i.e., unequal partitions of the cross sections and unequal heat transfer 
coefficients together with unequal capacity rate ratios of the fluid 
streams) but still for an infinitely large value of C,/C,;,. 

London has shown, as reported by Lambertson 1 , that a correction 


factor for conductivity in the direction of fluid fiow for approximately 


lNumber in brackets refers to bibliography on page 21 


equal fluid capacities, of the form: 


TT 11 


should result ín a somewhat pessimistic prediction of the reduction ín 
effectiveness. 

It is the objective of this thesis to develop and present a finite 
difference numerical analysis with solutions to solve for the conduction 
effect on effectiveness for the rotary regenerator. Since there are six 
dimensionless parameters required to specify behavior in the most general 
case, the ranges of the parameters covered in the solutions will be 
those of most interest to the gas turbine regenerator, however, the 


analysis will place no restrictions on these parameters. 





2. Method 

The conventional idealizations and boundary conditions assumed in 
the derivation of the governing differential equations, including the 
heat conduction iu the direction of fluid flow 4 , are the following: 
(See Fig. la) 

a. The thermal conductivity of the matrix material is zero in the 
direction of matrix metal flow. It has a finite value in the 
direction of fluid flow, and is infinite in the other direction 
normal to the fluid flow. 

b. The specific heats of the two fluids and matrix material are 
constant with temperature. 

c. No leakage of the fluids occurs either due to direct leakage or 
carry over, and each fluid flow is unmixed. 

а, Тһе convective conductance between the fluid and the matrix are 
constant with flow length. 

e. The fluids pass іп counter flow directions. 

f. Entering fluid temperatures are uniform over the flow inlet cross. 
section and constant with time. 

g. The matrix temperature gradient in the direction of fluid flow is 
zero at the ends, 

h. Regular periodic conditions are cstablished for all matrix 
elements, і.с., steady state condition. 

The differential equations derived by Schultz i21 are given in 

Appendix 1. 
It was suggested by Pevrberra amd carried out by Lambertson [1] for 


the case of no conduction, that a rotary regenerator could be divided 








Typical 


element 


FIG. la ILLUSTRATIVE MATRIX ARRANGEMENT 
AND FLUID FLOW 





into finite elenents as represented schematically in Fig. lb. This same 
method can be extended to the case with conduction, if the clementse are 
considered to be fixed ir space, cach element can be regarded as a cross 
flow heat exchanger with a gas stream and a metal stream. For an element 
on the side of Cmax», the heat transfer rate by convection is equal to the 


rate of change in enthalpy of the fluid across the element, 
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O а represents the mean temperature difference between the fluid 
and the matrix element. For a small enough elcment the arithmetic mean 


temperature difference may be assumed valid, so that 


αν = 2% +T = 3- 


4 \ 32 گے‎ x . = Ds k 


Г n=” 1 
| | 
| 
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u < 
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where the fluid and matrix temperatures are subscripted to indicate the 

average temperatures across the inlet and outlet of the clement for 

simple unmixed cross flow. - 
Considering an energy balance for the element, the energy transferred 

to the element by convection plus energy transferred to the element by 

conduction less the energy transferred out of the element by conduction 


must equal the energy stored in the element: 
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By using equations (1) and (2) the outlet temperatures for an 
clement on the sice of C,., may be solved for in terms of the remaining 
temperatures and the result expressed in dimensionless form (Appendix 1). 

Since the longitudinal heat conduction is zero at the ends, the 
energy balance equation is different from cquation (2) for the elements 
of the first and last row. Three different expressions for the matrix 
outlet temperatures are obtained: 


First row elements 


Middle row elenents 


مر 


Lsst row elements 


дмт Jj HL τι (5) 


By solving for the fluid outlet temperature on the side of Cray in terms 
of the matrix outlet temperature from equation (1), only one expression 


is obtained, which holds for all the elements on the side of C Rae! 


j F 7 a 4 / а. 


similarly, for av eicment on the side of Chin 


r 
{ fo ; гь 
E Е e — 7) 
{ : ou war u j 
= M a m М, i F ۲ | ( ғ 
| = 
with 
A (m) [ ۷ } = - 
we + 7 15 168: 
Ἢ, αἱ 1) Ж | А 
Energy balance; 
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\ Г H Е Ы 
— | e e т — — № Ш ы ید‎ E" 


y | 70 (8) 
By the same method the outlet temperatures of the elements on the side 
of Gain may be obtained using equations (7) and (8). 


First row elements 


I 
| | 3 i - τρ <. 


í Ч. д a ЛЕ = i i o ie (9) 


i У 
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Middle row elements 


i г. 4 |) | (10) 








Lost row elements 
} 
| - Р. (11) 


Solving for the fluid outlet temperature for the elements on the side 
of Cain from equation (7) gives 


Er: E | | En Eee | 
e A Lu, | & 


which holds for all the elements on the side Of Cmin. For the values 
of the constants see Appendix 1. 
In addition to the four dimensionless groupings for no conduction, 
Cnin/Cnax = capacity rate ratio of fluid flow streams 
Cr/Cmin = capacity rate ratio of rotor matrix to minimum 
fluid capacity rate 
ъл)“ = conductance ratio 
NTUo = . , over-nll number of transfer units; 
two additional dimensionless groupings are needed in the coefficients 
of equations (3), (4), (5), (6), (9), (10), (11), and (12) to take 


conduction into account. These conduction parameters may be defined as: 


\ — , conduction parameter on side of Cnax 


= , conduction parameter on side of Cin’ 


Also وی‎ 


since the thermal conductivity and the length of the matrix are the sace 


for both sides, and where 


i 
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- πετ» conduction area ratio. 


Another form of defining the conduction parameter ig 


πλ ہی‎ \ un ‚ total conduction parameter 


` یف‎ ` x 1 \, 


This form offers the acvantage that for a given value of the conduction 
parameter, the resulting effectiveness will not be affected by small 
changes їп As”. 

From the schematic representation shown in Fig. 1b, it should be 
noted that the left edge of the regenerator is physically the same as 
the right edge and, therefore, the matrix inlet temperature for the 
elements of the first column on the side of Cmax is the same as the 
outlet temperature of the corresponding element of the last column on 


the side of CURE Expressed mathematically, 


Sa 1 EF 3 №. +1) 


where i and f can take the values of l, 2, 3,............ i. 
This is referred to as the reversai condition. 

From equations (3), (4), (5), (9), (10), and (11) it is seen that 
in order to solve for the outlet matrix temperature for a particular 
element, it is necessary to know the matrix temperatures of the next 


element, An estimate of these values can be made by determining the 








temperature distribution for the case of zero conduction in tne direction 
of fluid flow, since the constants Dz aud Dg in equations (3), (4) and 
(5) and the constants F3 and Fp in equations (9), (10) and (11) are zero 
in this case. 


A temperature scale can be used for which the fluid entrance tempera- 


ture is zero at the side of C,;, and unity at the side of C... 


Ев (1,2) = 0 جات‎ lis dues «cilii e... «No 
tx(1,j) = 1 For | 1,2,395-........v. ا‎ 


For calculation purposes a temperature distribution is assumed on 
the left edge and the problem is solved for the no conduction case in 
order to obtain an initial estimate of matrix temperatures for working 
the problem with conduction. This also provides the no conduction 


effectiveness necessary for comparison. 


The outlet temperatures are calculated for every element by repetitive 


use of equations (3), (4), (5) ала (6) cepending on the Location of the 
element. The calculation is started with the first element, first 
column on the side of Cia, and then working down the column. When the 
first colum of clements is completed, the second column is calculated 
and so on until the outlet temperatures of all the elements on the side 
of Cmax are determined, It was found that convergence could be enhanced 
if two passes per column were made before proceeding to the next one. 
The reason for this ig that the equations for the outlet temperatures 
contain the matrix temperatures of the next element which are necessary 
to cstimatc. 


At the seal represented by the double line in Fig. lb, 








Tx(1,N, ^ 1) 7 Tn(£,1) 
since the matrix outlet temperature for the elements of the last column 
on the side of Cmax is physically the same as the inlet temperature of 
the corresponding element of the first colum on the side of Суд. 

The same method is applied to the side of Cmin 28 to the side of 
Cmax. If the temperature distribution assumed on the left edge vas 
correct it would then be duplicated on the right (i.e., the reversal 
condition is fulfilled). If, however, this is not the case, the 
resulting temperatures on the right side Тое н 41)! are now used on 
the left (T« (1,15) and the procedure is repeated. After each pass (i.e., 
the complete calculation of a temperature distribution) an energy 
balance is made; and before the solution is accepted for a particular 
get of parameters, the heat balance error together with the reversal 
condition has to be fulfilled to the specified accuracy. 


The effectiveness of a heat exchanger is defined (Appendis 2) as 
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From the conditions of the problem the above expression for effectiveness 


reduces to 
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The heat balance error was computed fron 
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That the solutiva will eonvergee (1..,, the cevepas! comditiour mt 
ami the ueat balance «rror converge to Zero) will depend on tae number 

of elements used, More important than this is the relation existing 
between the nuab ir of subdivisions for the three streams. The sufficient 


but not necessary conditions for convergence are (Appendix 1) that 
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Fig. 2 is a plot of the ditferent values of effectiveneag obtained 


for the following set of parameters: 
Cuial Ema: * 1-03 Cr/Guin 7 2.0, (ВА)* = 1.0, NTUo = 10.0 \=0.1, As = 1.0. 


Seviral curves are obtained depending on the relation existing between 
the aumber of subdivisions for every streem, To extrapolate the 
effectiveness to an infinite number of elements, these curves may be 
approximated by straight Lines. The error introduced will end on 
the nurber of elements used and will not be constant for all the range 
of the porzneters. As an illustration consider the sot of parameters 
used to plot Fig. 2. Three cases are investigated: 

1) Y, Ҡа Мұ/2 

2) Me = № = Ne 

3) ty =u, » 2N, 

A lincar extrapolation is performed using the values of effectiveness 


obtained for values of 1 - , From the results it is found 


Nr (Nyx + №) 





that the minimum Ny; necessary to give an error in the extropolated 


effectiveness of less than 0.01% is 18, for which е. = 76,577 
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for the three causes. The determining factor in the extrapolation is 

the number of matrix elements. The higher the value of Ny used for the 
extrapolation the smaller the error in the straight line approximation. 
Unfortunately, an increase in Ng will bring a higher increase of N, and 


Nn, Since the convergence condition for this example is 


That this is only a sufficient condition can be seen from the case of 

Ny = Ny = Np/2 
for which the solution converges for values of Ny >5; and when 

Ny = N. = № 
the solution will still converge for Ny, >20. It is true, however, that 
the closer the number of matrix elements are to the sufficient convergence 
condition, the fewer passes it takes for the solution to converge to 
the specified error in the effectiveness. Table 1 shows the number 
of passes necessary to obtain a .olution with a heat balance error of 
less than 0.005%. 

For a given set of parameters the effectiveness was determined 

first using one set of subdivisions and then recalculated increasing 
the number of subdivisions. 4A linear e:xirapolation was performed to obtain 
the effectiveness corresponding to an infinite number of elements. The 
accuracy for each value of effectiveness was spccificd to be of four 
significant figures. When these values were extrapolated to an infinite 
number of elements, the accuracy was reduced to only three significant 
figures since a compromise had to be made with the time available for 


computations and the number of runs to ос made. However, the nethol 
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Number of passes necessary to obtain an effectiveness with a heat 


balance error of ices than 0.C05 for the 





following set of parameters: 








ЕХ: Е10,Сс/01, 52.59, (М) = 1.0, 
NTUo = 10.0 , \= 0.1, As” = 1.0 

- pem | 
Ny Ny Nn Mo. of passes | г ей 
12 12 12 27 76.581 
B? 25 24 (2 76.581 
16 7 7 41 
14 14 14 34, 
14 28 28 27 
16 16 76.574 
16 32 76. 574 
16 
13 18 
10 36 
20 2 76.572 

Ec 22 2 71 76.571 


* о„ is obtained by a linear extrapolation of the effectiveness 


for Ny and N, -1. 





imposes ro restrict,on ак to the ultimate accuracy of the solution. 


When working with very high valucs of CARR tne computation time 
becomes considerably long. This is due to the fact that the finel 
temperatures differ very little from the assumed initial temperatures and 
a large number of elcments is required. The starting estimates of matrix 
temperatures then become very important. For the Limiting case of C,/C,i, 
equal to infinite, there is no temperature variation in the direction 
of metal flow so that the method without modifications is not applicable. 
For this case, however, the exact solution of Hahnemann EX is available. 
The Fortran system was used to program the above method for the 
CDC 1604 digital computer. The matrix notation employed in the equations 
derived previously is very convenient in writing up the program, since 


several loops can be made by just changing the subscripts of the variables 


involved. The complete cowputer program is given in Appendix 2. 
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3. Discussion of results 

The solutions generated are presented in tabular form in Tables 2 
through 25 and graphically in Figs. 3 and 4. The accuracy in the 
effectiveness was specified to be of three significant figures. It was 
found that for 

0.7 (hA)*< As* <1.3 (hay 
the resulting values of aci vere almost ídentical to the ones 
calculated for 
As” = (hA)" 
so that these were the values used in preparing the Tables. 

By using NfUo and Л аз defined in the nomenclature, which contain 
(ҺА)* and As” explicitly, the influence of these two parameters on the 
conduction effect, =>, were founá to be very small for values of 
Cnin/Cmax > 0.9 and 130.1, which is the range for the gas turbine 
regenerator problem. 

Figs. 3 and 4 show the effect of conduction on effectiveness for 
values Of Cmin/Cmax = 1.0 and 0.9 respectively, with different values of 
NTUO, Cr/Cmin and \, up to a limit of 10%. The conduction effect 
continues to increase with increasing A. For small values of A this 
effect increases nearly exponentiolly. By further increasing A, a 
saturation is reached which approaches the limiting condition of \ = оо. 
Fig. 5 is a representative curve which shows the limiting influence on 
the conduction effect by continuously increasing A for a given set of 
parameters. 

The conduction effect for values of CL/Oq,4 = 10.0 are almost identical 


with the values reported by Schultz [2 | and Hahnemann [3 | ; [6 | for the 


imiting case of C./C,;, = 09. 


م 
c‏ 








For small values of NTUo the effect of increasing Cr/Cain 13 to 
increase slightly the conduction effect, as can be seen from Fig. 6. 
For higher valucs of NTUo the effect of conduction is reduced with 
increesiug CCA > 

The effect of increasing NTUo upon ac is intimately related to the 
value of C-/Cmin 20d Cnin“Cmax: This effect can best be explained by 
observing Fig. 6. For the case where Cmiín?Cmax = 1.0 the effect of 
increasing NTUo will result in an increase in € regardless of the value 
of C./C in; Of course, the higher the value of Cy,/Cgin the smaller this 
effect will be for relatively high values of NTUo. Now, by reducing 
Cmin/Snax and increasing C,/C,;, a point is reached in which an increase 


in the NTUo will result in a reduction in =. The point where this 


с 


occurs will depend upon the value of the conduction parameter A. 
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а. Sumsarv end Conclusiona 

A finite difference numerical analysis {s presented which will 
Getermine the effect of the Longitudinal heat conduction in rotary 
regenerators for steady state conditions. 

The ranges of governing parameters which have been covered are: 


0.5 «C C LO 


in?! шах © 
1.0 Str боза 210.0 


1.0 < i#fUo < 20.9 


0.15 € 1.3 


The ranges of parameters of interest for gas turbine regenerator 
application are plotted in Figs. 3 and 4. A more accurate determination 
of the conduction effect than that given in the figures can be made by 
interpolating values from the Tables. 

The computer program used to carry out the calculations on the 
CDC 1604 dipital computer is included in Appendix 2. This program 
imposes no restrictions on the parameters for which the analysis would be 
applicable. It is recommended, however, that for values of Cr/Cmin >10 
the solution for C./Cmin = 10 be used. For this case the computer time 
for a solution increases considerably and the error introduced by assuming 


Cr /Coin = 10 is negligible. 
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Ь 7, i 
NTUo | o 0 
ΠΠ | ο 
2 | 60.06 | 0 
3 | 66.72 0 
4 x 70.86 | 0 
| 5 73.75 x 0 
6 175.92 ' 
7 | 77.63 | 
8| 79.01 | 
9 80.16 x 
10 |81.15 |1 
B 82.74 
14 86.00 
x 16 | 85.01 
18 | 85.86 
20 | 86.58 | 1 
x 40 |90.37 | 1 
| 100 93.87 |2 
| 500 : 97.87 | 4 
| 
x 1000 98.81 x 4 
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891 0,90 1.79 3.49 6.59 11.74 ob 
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* 
Onin/Onax 7 .9 , Cy/Gain = 10.0 , (hA) = А8 


سے ےہ ےہ ہے کے ےس تہ “τ᾽.‏ ہو ا 


Conduction effect —— % 


------------.------.. әмме > q --- 


وک ہے ہے πι.‏ 
NTUO \‏ 





TABLE 14 





86. 


a4, 





—_ —n —— 


. 01 
0.70 1.39 
0.75 1.48 
0.65 1.33 
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. 04 


.45 
. 06 
. 38 
«96 


‚68 


. 88 


‚70 


m U = 


.08 


2.56 
3.77 
4.39 
4.76 


5.00 


5.48 


5.35 


LONGITUD, ERAT CONDUCTION EFFECT 


4.23 
6.49 
7.68 
8.39 


8.87 


9.87 


10.03 


A s s Que me mcm 
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| TABLE 15  LOSGLITUDINAL HEAT CONDUCTION EFFECT 
| 
* 
| Саіп/ вах = -9 » Cr/Cain = 1.0 , (hA) = As“ = .5 
Conduction effect کے‎ 6 
2:٣ ө | 
.03 .06 2 . 24 .48 
„03 1.88 3.23 ӘСІП 7.23 
43 2.67 4.76 7.89 11.81 
69 3.17 5.70 9.55 14.48 
„88 3.55 6.40 10.73 16.29 
5 75147 1.06 2.05 3.86 <.94 11.63 17.64 
10 | 82.98 | 1.35 2.61 &.89 8.69 14.33 21.38 | 
20 88. 34 1.68 3.23 5.97 10.39 16.70 24.35 
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ТАЛЫ. 16 LONGITUDINAL 


| 
| 
Cain’ max ای‎ Cr/Cqin 7 2.0 , (hA)* = As" = „5 
| | mo Conduction effect, 26 "s 
N | 0 03 12 21 30 
i120 ° : se 4 
<= | AN A о وھ‎ 
1 | 50.27 % E 3.40 4.96 5.09 
E 2 x 66.88 | 1.53 5.08 7.71 9.73 
Е | 75.23 | INE 6.01 9.23 11.76 
Ж. 20.32 | 1.93 6.63 10.24 13.08 
| 
x 5 83.77 | 2.05 7.09 10.97 14.03 
6 | 86.28 x 2.14 7.45 11.53 14.75 
l 
7 x 88.18 | 2% 77 11.98 15.32 
3 3 89.69 x 2.29 7.97 12.34 15.78 
9 90.90 | 2.34 8.17 12.64 16.17 
| 
x 10 x 91.91 2.38 8.34 12.90 16.49 
12 x 93.47 2.45 8.60 13.31 17.00 
" | 94.62 2.50 8.80 13.61 17.39 
| 
| | 
| 16 | 95.50 2.54 8.96 13.85 17.67 
| 
18 | 96.20 2.56 9.07 14.03 17.90 
20 | 96.76 | 2T 9.16 14.18 18.09 


< 
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nn nn ees t Y=" —‏ - ریو ےم — س 


TADLU 17  LOWGI'TUDINAL AZAY COUBMCTEGR ErrscT 


«ᾱ- 
а еъ 15% ve (u^ y ہے‎ 
: "s t = / = ^ = 
ап a E mt 7 139.0 , (Δ Ав . 


о a د کے چ چو‎ Me s чан 7 En PS SY SG AS ےھ لے‎ 


| | еы Y 
8% | Conduction effect  —— % 
| | 


0 | 0.005 0.02 0.06 0.12 











0.80 Pw 2g ἌΠΟ 


0.92 1776 3489 6.05 


1.08 E299 0ء8‎ 6.39 


ہے ا ہے ہے ص ہہ < بے arene ee ee M et wm‏ م د مہا ا س A ee a ee‏ 


| 
| 

82.98 0.98 1.91 3.62 6.58 
| 
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=. ema meinamme A s «do amend douce t are aom m e بج وت‎ 
π.μ μμ μμ og HP I AN. O UR B pln 60 7 


9, ИА "1.0, (зА)" = а= .25 
Conduction effect Р == ۸ 

0.025 0.05 0.1 0.2 9.6 0,8 

0.02 1.62 2.90 4.62 TOI 9.57 

1.18 2.26 4.13 7.08 12.98 15.09 

1.39 2.65 4.87 2.29 13.10 18.12 

1.54 2.94 5.41 9.32 14.54 20.11 

1.66 3.17 5.83 10.03 15.61 21.54 

2.05 3.92 7.0 12.16 18.61 25.35 

2.43 4.67 0.45 14.04 21.04 22,16 | 
ا ا ا‎ m nn nn nn ШЕ. | 


— س ےس‎ Ыыы. o ص‎ MM جد کے ب ااه‎ So mee а н а nm س٢ہ‎ 





m nome سے‎ m am EU om سے‎ e LE e — - 


„ VP OPCTIBDIVAL "EAT UoUDUTTIOS ZDPPECT 


E е | е 
изу ο, mney 7 10.0, (М) -4α = «25 


4 $ ΄ 


ΗΓ a —‏ دم  _ a.‏ س 
A‏ 
A‏ 
مہہ 


To‏ —= مہ یت یو رٹ 


meee жыз À— σσ. -- — —a e. . peat سس‎ ——— — «ρα 














v. я 0.05 0.1 072 0.4 6.5 


«τρ‏ حسم ھہ 
e en 9À -. = ——————————————— а 05 — = P AIMO. rfe Rss o A en ran‏ 09 ہے — À—‏ ےہ ہے MÀ‏ — — 


0 43 L.7é 3.19 5,1» 3 10. М 
1.23 2.45 4.49 7.72 1290€ 16,74 
1.47 28? 1.21 9.06 14.32 20,123 
36 3.03 7.64 3.87 Lon есе 


۰ھ 16.62 10,41 5.93 3:17 1.64 


l./3 3.41 6.28 11.56 18.67 26,62 





29 3526 6.76 11 ag 19.54 28.00 
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= а= .. — — — ap t m et m m a — 0 س س ا‎ 


TABLE 20  LONGIIUDIMAL HEAT GONDUCTIOR EFFECT 


س س ی ا 


% ` , — ہما 














E ex ا‎ Conduction effect, = 7 
0.01 0.02 0.04 0.65 0.16 0.32 
0.35 0.68 1.29 2.37 4.15 στι 
0.41 0.80 1.53 2.84 5.02 8.25 
0.45 0.89 1.71 3.19 5.65 9.29 
0.49 0.97 1.36 3.46 6.12 10.06 
0.55 1.26 2.39 4.38 7.57 1215 
0.86 С 3.04 5.37 8.94 13.97 | 











—— αν со 
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-> u un i = - - - ev + з». --- — — — ووس س‎ = — — ------ - e ae м 


тазы Z1 LONGITUDINAL HEAT CONDUETLON EFFECT 


——— 











0.01 0.02 0.04 0.08 022 0.32 
0.25 0:53 0.98 1.74 2.90 4.37 
0.37 0.71 1.35 2.49 4.36 7.06 
0.33 0.74 1.43 2.68 4.82 9.08 
0.35 0.70 1836 2.60 4.80 8.29 
0.21 9.62 222 2229 4.52 9.1 

0.10 0522 .49 1.14 2.68 5.96 
0.01 0.01 0.05 0.23 1:02 3.65 
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DO ome dp ma 








| Gin qim 0.5 , Co/Ggn S O, aaa ОБ 
x Е% | Conduction effect, &— Za 
x | 
TIO EX ‚ 9.01 0.02 0,0% 9.03 0,16 0.32 
Yo ог. ‚53 0.98 172 2.75 3.99 
| | 
24 065.501 0.35 . 69 [130 2,40 жш © 6,50 
359 72.44 Ü оз ‚78 1,49 2.78 4 ge 8.04 
| 
^ | 76.43: 9.52 ‚87 ΠΠ 3.04 5.44 9.04 
| | 
S 79.10 | 0.45 0.89 1.72 3.25 5.86 9.79 
| 
10 | 85.44 | 0.56 1.09 2.12 4.09 78 11.95 
| | 
20 | 89.72 | 0,7% 1.54 2.78 4.95 8.58 13.90 
| | | 
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m‏ ے۔ -ہو بھی جم 





%9 


ee کے ص سے‎ ае das ee TA A AA PILA AU түзы» IA A =... - 


[ 
| 
| 
| 
| 


MÀ m m ου  —Ñ.—OÁo D ———.— sn —vvnI O — A I 2 


AAA ARIES ns ea ee ο 


—— — 5 — Le SY STS — — s... -... v. — ےچہہ‎ ыт = 


۹ dita ut 
x) 4 т, “4, 1 + δὲ ο: {5 "uc Ж ОҢ 
y 
= ۴ھ‎ F e 
"94 5 3l + $ (Е) мы Oe E 


a 





Cunduction etfe E gF 7 


"` a —... nn A πι πι η 





№. < 0.04 0.98 0. 5 0:22 


;. 61 ιά 1,95 3.1% 4.53 
2.33 1.27 7 fj m 7.76 
0.37 1.66 3.1. 5.61 9.12 
0,12 1.62 3.10 m 9.75 


9,74 1,47 2.90 4.02 0.76 


ο 0.83 1.54 3.74 0.13 


* 
——.—.— سے سس سے وم سب ہے 
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TABLE 24 LONGITUDINAL PEAT COUDUCYIGN EFFECT 
е % | Conduction effect, — А 
=| u A mx 
0 | ‚ol 02 ‚04 OR „16 ‚ 32 
50,21 | 6.34 0.65 1.13 2.01 3.12 e 
54.39 | 0.38 0.75 1,42 2.61 4.47 6.89 
70,92 | 0.39 0.77 τ" 2.89 5.19 8.88 
{ 

74.91 | 0,39 0.79 1.58 3.09 5.70 9.43 | 
77.57 | 0.40 0.81 1.64 3.26 6.11 10.24 
94.12 ! 0.49 0.99 1.99 3.97 7458 12,71 


88.72 0.56 1,39 2.54 &. 92 9.07 14,97 
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TABLE 25 


a > 


= 


.- — ہے — 


LONGITUDINAL HEAT CUsDUCITION EFFECT 


PI — x. nn ULA nn озше»‏ ےم ہے 


æ smm — „4 


д ” 
и бак 05, С/С = 10.0, GAD hS Ls | 


.02 


0.81 
1.11 
1,18 
1,12 


1.02 


0.41 


0.04 
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Conduction effect, = 


‚ 04 
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2279 
4.24 
4.28 


4.10 


2.59 
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6 „32 
3.92 5,56 
6.26 9,59 
282 11.66 
7.77 12.62 
7273 12.99 
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4.29 10.59 
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APPENDIX 1. 


herivation of iterative equations.- 


Equations (1) and (2) from page S can be put in the following form: 


ПЕ. ® С | 7 a π A‏ تی 
Р w | E i l ux : `‏ 
“ер” thue ο T uen ap T "egy еде) |‏ 
B € » ^1 ADN. t (3 Ius стл‏ 
Ce m f ат‏ — 
νο "m с j ) | .. Ma: 1) ST 4 )‏ 
E А‏ 
Pas = \ — i‏ 
ДЕ (А T sp | т‏ == 
A ) 541: ) X UT 1) Те J‏ 





where 
N i WES | 
Ке ae ای‎ 
А. ` PI 
-- 
К. کے‎ — = == joa 2 ےصح‎ = 
and 
а = \ A ee | | 7 ‚ 
Equations (1.1) and (1.2) may be written as 
\ P x 
M e c τ. 
Ke МЕ ` 
and 
NL к ۱ и 7 NA Е. A 
سر‎ | | ΤΕ L 
` \ š n NIQ τ `4 
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(1.1) 


(1.2) 


(1.3) 


(1.4) 





ру adding to equation (1.3) the product of (1+C, ) (1.4) and precia ef 


4 
| | (4) 
where 
and A c ` 
From equatíou (1.1) 
| o Γ ! (6) 
P 1) 
where 
٦ 
/ e and = 
For the elements from the first and last row equation (1.2) becomes 
[ " i I hr 1 5) 
and 
h 1 2 (1.6) 





resoeclivel,. I. t.e same procedure as before and solving for the matrix 


outlet temperatures, the following equations are obtained; 


and 


for the first and last row elements respectively, where 


LE 


` м and |) a. 3 


Equations (7) and (8) may be written as 


and 


Equation (1.8) reduces to 


л 
ғә 


(3) 


(5) 


(1249 


1 5 


- == 





A 


f um Wk E 
d = Ms 
TN- 4) AN + [ye ) R| btu, asi Για М ч 


апа 


/ (1.10) 


By comparison of equations (1.7) and (1,8) with equations (1.1) and 
(1.2) it follows that the form 15 the same and only the coefficients are 
different, The same holds true for equations (1.9) and (1.10) wíth 


equations (1.6) and (1.5) respectively. By substitution the following 


equations are obtained: 


۰ یش سس 
A 1 |‏ 
Г‏ 
UN к (аі a» '‏ 
Ley -‏ - 


(10) 


(11) 


Ci 
ч 


for the first, middle and last row elements respectively and 
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(14) 


which is the same for ail the elements, where 


ہہ — 
— 


атат 


and 
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COMTCE RENCE Condition,“ 


Equation (1,2) can be put in che relloving form: 


= B. (1.11) 


n tha energy balance equations, the side of C425, was assumed to be the 


hot side go that by the second law of thermodynamics 


and from equation (1.11) 


(1.12) 


ہے 


Reaxrii n^: tos equation (1.12) 


(1.13) 


Again by the second law oi thermodynamics the matrix temperature variation 


іп the direction of fluid flow is such that 


| (1.14) 





апа 


(1.15) 
So that 
(1.16) 
Substituting equation (1.16) in equation (1.13) 
۹ — — | (1.17) 


Three cases are investigated: 


l.- TÉ с, с, then equation (1.17) becomes 


ФУ I£ Ca = O then equation (1.17) becomes 


v 14 


3.- τε с, - 201 = О then equation (1.17) becomes 


From the above cases it can be seen that a sufficient condition for 


equation (2) to converge is that 


Substituting the values for the constants and rearranging yields 
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For the first and last row elements the convecgence criterion 


becanes 


which for substitution of the parameters gives 





pA 
A 


” 


This condition is less stringent than for the middle row elements. 
Similarly, for the side o. “>in the sufficient condition 


for the middle row elements becomes 


which for substitution of the parameters gives 


For the first and last row elements it is found that 


which requires that 


is the sufficient condition. This again is a less stringent condition 


than that imposed by the middle row elements. 

















2iifeTential equaoiio ü.- 

Hausen 4 has given д theory for regercrátor$ without neut 
conduction. When the same nssureptions ¿re made, but the longitudinal 
conduction of heat is included, (See page + ) the differential 
equations obtained are 


for the side of € 


max (838 flowing from . = 0 to = L); 


and for the side of Gain "reversed gas flow) 


A temperature scale can be used for which the gas entrance 
temperature at the side o! Cain is zero and that at the side of C 


вах 


is unity. The boundary cot.J.tions are then 


ےمذ —— 
سسسوے سے — 


(9 


u. 





APPENDIX 2 
Computer program. - 

The Fortran system vas used to program the problem for the CDC 160% 
digital computer. 

The sis: dimensionless parameters are used ag input together with 
increments and factors to change the values of NTUo, and кй 30 that 
solutions can be generated for several combinations of the parameters 
without having to feed into the computer these different seta of para- 
meters every time. 

Another input is the number of subdivisions of the three streams 
together with the increments, so that when the computer finishes the 
first series of runs for a given number of subdivisions it automatically 
increases the subdivisions by the specified increments and the two values 
of effectiveness obtained are extrapolated to an infinite number of 
subdivisions. 

In order to obtain the initial estimates of matrix temperatures, 
the problem is first solved neglecting the longitudinal heat conduction 
as explained previously. For this case the equations derived in Appencix 1 


are simplified to (for detailed derivation see Reference 1) 
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MM کم‎ 
a> 
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ZI í 

Ú и 

гү 
5 





The inlet watrix tesperatures necessary to initiate the 
problem for zero heat conduction are estimated by assuming |, 5 
(See Tig.1:) and feeding this value into the computer. The remaining 


inlet matrix temperatures are obtained from the relation 
= | — 


The number of passes necessary to meet the desired accuracy in the 
solution changes accordingly with the initíal estimates of matrix inlet 
temperatures, but considering that the computer time required to vork the 
problem for zero heat conduction is so small es compared to the case with 
conduction, that these estimates do not become important. Only for high 
values ot C; Cuin these estimates may be of significance since for this 
case the change in the reversal condition after cvery pass is very Email. 

in general, it was found that an estimate: T, (1,1) 


of 0.4 to 0.6 £or C-/Cmin 3.0 


and 0.7 to 0.9 for Cr/Cmin >5.0 


together with the above relation for T4441) should converge to the 
solution in a reasonable time. 

Tne temperature distribution is calculated by repetitive use of the 
iteration equations. 


The effectiveness may be defined as 


Бас: "πι 
Or in oe — Ui) Ll. lna 
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t9 CAE conditions ol δὶ negok lai 


| 


1 

| 
SA 

ñ 





ہے 
z|‏ 


therefore, the expression for effectiveness for the definition reduces to 


М, 
A > ia i 
n — 


vhich is the expression used for the computation of the cffectiveness. 


The heat balance error is computed from the relation 





A ` حا‎ Te 22 t C ۲ 
error = ee پا ہب‎ πι 
í f = \ 
EUNT ` A ) 
which reduces to 
error = > سے‎ 
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Tac extrapolated effectiveness is computed fron che relation 


š rs e 
(«зс : de . 0 


«леге 


КЕ فو راج‎ A. AA 


The input parameters in Fortran notation are 


Pl = Cain Овах › F2 = С/С . 23 = (A) , 


p4 = NTUo , PS = \ , рб = Аз 
The initial values of NTVo, X and As” are changec by the following 
relations 
P4 = P4 + PALIN 
P5 = P54P5FA + PSIN 
p6 = p6 - Ρ6ΙΝ 
where PAIN, P5SFA, SIN, POIN are also input values. The inputs P4FI, 
PSFT and P6FI are the final values of P4, PS, and 56 respectively. 
The program is written so that the computer after working the 
initial set of parameters, i’ ¿c 2. cos ` by P6IN every time until 
рб becomes less than P6FI. Then it sets F6 equal to the initial value 
and changes PS. The procedure is again repeated until P5 becomes greater 
than P5FI. Then P5 is set equal to the initial value and P4 is increased 
by P4IN. The cycle is repeated until P4 becomes equal to or greater than 


FPL. 
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= 4 pt $ + ہہ‎ 
Ж. `... ж - - т = ۶ ‘ 
burlny Che гет ri La "м ССА suc madé by Phe Ця ot Pi 


e 
= 


. - X I" ex ہے‎ 1 1. = е” ج‎ 6 1» P . t ` 
СВЕ нет са сіз s ul rz со walk oF ach 24cis.2no іл usa Y a Leas 


паре or ofS mo convorge metan tle maxima nome Ql ow npa ؟۶‎ 
Жа, a i - ге. p x He gs г, Аға eh 
5 pP: ۹ Y 1 + [i * 8 4 ' “ “ 78 є? 461 $ ha Show Ene compl c 1 E 1. 201 T e 1 vt & = м» % 11 


κ. * fe Ἂν i ae | > А 4 سے‎ > j . ь ہہ‎ - 
¿ogro which le folloucd By the Fontan listing of ch Brom. 





| 1 INPUT ШІН 


| M = PASS NUMBER 


e COMPUTE INITIAL Tx(1,1) | IC = 1, INITIAL SUBDIVISIONS 


| 1С 2, INCAE 'X-NTED SUBDIVISIONS 
| DEESSET tx(1,j) = 1,0 AND tn(1,8) = О 


NO CONDUCTION CASE 


st 













COMPUTE Tx(1,3) , ЩЕ к) 


CONDUCTION CASE 


| 23 COMPUTE CONSTANTS 
| 24 SET M = 1 

25 ITERATION ON SIDE OF Cray 
26 ITERATICN ON SIDE OF Cmin 
27 COMPUTE € AND ERROR 


» € AND ERROR 











YES ест جح‎ 








& M£SPECIFIED MAXIMU 






YES 
(€ - €, _, ^ SPECIFIED m | 


ہیں 















35 CHECK CONVERGENCE 
AND PRINT 


IC 


" 

P 
m 
O 
u 
гә 





36 INCREASE M BY 1 


38 REVERSAL CONDITION 






3i COMPUTE OG. 
AND PRINT 


30 STORE € AND ERROR 





37 SENSE SAITCH 





REPLACE Tx(i,1) BY 







Tn(f ,Nn+1) 





|39 PRINT € AND ERROR 






32 SENSE SWITCH 


LO SENSE SWITCH 







33 PRINT REVERSAL 
CONDITION 





41 REPLACE Tx(i,3) AND Tn(f,g) 









BY INITIAL ESTIMATES 





+ 
42 REDUCE As 





NO YES 


As < FINAL VALUE 












р 
144 SET as” = INITIAL VALUE | 
1 


INCREASE Aa 


— o a 


| 


YES S 








| 
| 


[52 sense sarren — 6) 


53 INCREASE 






48 INCREASE 









5) SET INITIAL 





Nr. Ny, Na 


(DH so ser ιο 


Np Ny. Np 





№, Nx, Nn 





54 SET IC = 1 


FIG. ?a GENERAL FLOW DIAGRAM 


о 
Qi 





4 COMPUTE TANSTANTS 
S SET É < 1 
"6 ITERATO: a SE oF Coax 
| 7 ITERATION ON SINE OF Cain 
|8 COMPUTE @ AND ERROR 





ERROR $ Si TCIFIED AXIOM a = 


= 


——  - — ο. 


T" TA E YES " Ж 
NO m 
С re ,) € SPECIFIED MAXIZUR 
м M- 
TE = 1 






| 11 STORE 











| COMPUTE е, 15 CHECK CONVERGENCE INCREASE M BY 1 
| i AND PRINT AN INT 
Tx(i,1),€, EROR ' Ξ E 


| DLE σα 


18 REVERSAL CONDITION 


|. (52-117 әке switch | o- "RPLACE DECEM) 


BY Tn(f,Nn+1) 
| ;——— e — J 





i 


К PRINT € AND ERROR | 
| 


| | 


(25) (5)—{ sense SAITCH 0 


| 


— r  —ns——OF? سے‎ 


22 INCREASE Np, Nx, Nm 


FIS. 20 FLOW DIA1«AM FO3 NO CONDUCTION 
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EE 





PROORAM REMARKS 





PROGRAM BAHNKE — 
OIMENSION TH(bQ, hO)UH(& O, UOT, TC(Lu0,uUOT,UC(S 0,5 0) 4 ERR(9 2) ,TCOqO92), 
IUHC( 405 401 ,UCCCSO0, hQ) UL T(&01 ,UH2(80) ,EFCÉ20,20) ,ERT 20,201 
FORMAT(6F5.21 
REAO 1,P1,P2,P3,Ph,PulN, Pufl 
2 FORMAT (6151 
READ 2,NM,NH,NC,NM],NH(SNCI 
FORMAT(FS.3) — 
READ 3,UH( 1,11 
ЕМеММ 
00 401 K=1,NM 
ЕЕК 
UH(K,11=((EM+1.0-E le UHI1+1)1/EM 
پا‎ FORMAT(7F 7.21 
REAO &,P5,PSIN,POFA,PSET,PO,POIN,PÓFI 
S5=p5 
S6=P6 
NM2=NM+ NM] 
NH2=NH+NHI 
NC2-NC*NCI 
NE I NM» (NH*NCI 
RE 1=NE f | 
NE2zNM2* (NH2*NC2) 
RE2-NE2 


= E С 
=R - 1 














INPUT 








NM ж Nr, NH = Nx, NC я Ид 


















INPUT: UE(1,1) » Tx(1,1) 





ESTIMATES OF Tr(1,1) 


STORE INITIAL P5 AND P6 


NUMBER OF SUBDIVISIONS FOR end RUN 


























RE = VARIABLE USED TO EXTRAPOLATE EFFECTIVENESS 


EM2-NM2 

LOS L=1,NM2 
tst 
RT=1.0-((E-1.0)/EM2) 

ὁ UH2(LI=UH( 1, 11oRT 
FORMATL//IUHNO OF ELEMENTS3X,3IS UK, 305) 
WRITE QUTPUT TAPE 2,/,NM,NH,NC,NM2,NH2,NC2 
PRINT. 7 ,NM,NH,NC ,NM2,NH2 ,NC2 


8 1C=1 INDICATES let RUN, USING INITIAL Nọ OF SUBDIV, 


FIXED TO FLOATINO POINT CONVERSION 


ESTIMATES OF Tz(4,1) FOR 2nd ВОИ 











NMMzNM-( 
NHPzNH*( 
NCP=NC+1 
NMP=NM+ 1 
00 10 J=1,NH 
10 TH(1,J1=1.0 
DO 1 






DEFINITION OF NEW VARIABLES 




















SET tx(1,3) a l.O AND tn(1,g) = 0.0 


12  ATO=(.0+EM/(Pler20EH1+(2.00£MeP31/(PuePle().0+P31) 
A1=2.0/A10 
А20<%.0%(Р1еР2еЕНІ/ЕМ»(2.0еЕнер3ер2і/(Рие(1.0%03511 
A2=2.0/A20 
A30=1.Q+EM/(P20EC1 +(2.00EM)/(Pue(1.0+P311 
A322.0/A30 
ALO=1.0+(P29 EC 1/EM +12.09ECeP21/(Puo().0+P31) 
А8-2.0/А%0 
с HOT SIOE NO CONDUCT) ON 
121 Mal | M DENOTES PASS NUMBER 
13 00 )5 J=1,¿NH 
UO 15 1=1,NM 
THUS GJI)STHOL, JISATe CUL CL ,JI- THO S411 
کا‎  UHLI, Je l)sUHCLE, J1-A29 CUR CL ,J1- THCT JD) 
COLO SIDE NO CONDUCTION 
DO 17 к=1, NM 
UC(K , 1H23UHCNH& TI-K 1۱ء‎ 
DO 19 L=1,NC 
0019 K=(,NM 
ТОКИ, Е) =ТС (К, AS (UC UK L2 -TCEK SL 1) 
UCLK,L* 1H-UC(C LL I- Ab CUC CK  L) Το (κεί ) 
G1z0.0 
00 20 J=1,NH 

20 СТЕТНОММ+ 1, 1+5 1 
THO=G1/tH 


CONSTANTS FOR THE NO CONDUCTION CASE 


ITERATION ON SIDE OF Coax 










ITERATION ON SIDE OF Cain 


СТЕГІ 


оү 





o 
1 








mu τει 


1 








G2=0.0 
Do 21 L=1,NC 


----------- 











21 G2=TC(NM+1,L 1+G2 
_  __TCO(M)=G2/EC 
ERR(MI=1,0-(P1eTCO1MI1)/(1.0- 1110) 
IFCERROM) )25,24,2% 
23 — ERRIMIS-ERR(M) 
24 — IE(ERRIM)-5.0F-61 30, 30,25 
| 28 (Ғ(90-М1 206,26,215 I 
| 26 IFCERR(M-1)-ERRIM)) 29,29,27 q . 
| 27 PAUSE 11 
| 28  IFISENSE SWITCH 11 32,121 
| 29 PAUSE I2 
ا ا س‎ E 
' 215 DO 22 I=1,Nm 
22  UH(I+11=UC(NM+ I-1,NC+1) 
M=M+ 1 
60 το 13 IAS 
30  T2=TCO(M-I)-TCO(M) 
IF (12) 301,302, 302 
$01 τ2--τ2 
302  IFUF2-4.0€-61 305,305,25 
305 1FU(C-1) 31,31, 33 и 
31 LO 312 J=l,NM 
312. UHILJ)SUHCJ, 11 
EFL, 1)=TCO(MI 
ER( 1, ()=ERR(M) 
"7 07ص99‎ КО зм... 
32 ҒОНМАТ(%ҒВ.2,(5,Ғ10.6,Ғ9.2,151 


РНІМТ 52,РІ,Р2,Р3,Ри,М, ТСО(У?,ЕЯК(М?, (С 
PAUSE 13 
IF(SENSE SWITCH 


321,121 


1 





| 321  IF(SENSE SWITCH 2) 1,322 0 
322 NM=NM2 
| NH=NH2 
| NC=NC2 
loros. "D ه۷"‎ ma. E 
33 — EFFeTCOLM)e(1.O*RE)-EFCT, T len 
عند‎ ی٣ج‎ --- - 
Із  ҒОЯМАТ(//%Ғ8.2.16х,21Ғ10.6)9.21,Ғ10261 کیہ‎ 


| 


| COMPUTE HEAT BALANCE ERROR AND COMPARE IT 


‚ WITH MAXIMUM SPECIFIED 


J- - 








WRITE OUTPUT TAPE A SEI O аен, тураа 


00 35 L=1,NM 7 7 


35 — UH2(LI=UHIL+ 1) 
E 1712 2 ( 
δι .. 

DO 44 I=1,NM 
00 42 J=1,NHP 
UHC(I+JISUH( TL, 9) 

43 00 && L=I,NCP 

uu UCCU(+LISUCIIALA 
ς CONSTANTS 

50 Еітіріе(1.0%Р322/(2.0еЕУ) 

1٤ ۶١ ء١۱‎ 
F2=(P2eECI/EM 
٤320.5105 
EU=El+(ENIe(2.09E3+£21) 
E5st Ie(ETIo(E2«t 3) 1 
F1s02.0«E11/Fh 
F2=(El+(Elle(2.0*E3-E2)))/£u 
F3=(E3*eE 11)/Eu 
Fus(1.0-ETVI/E 11 
۲۱ 
۲۵۰۱20 ۰۹ 
FT=(EleCElle(E3-E2J11/ES 
Εβ-{Ε5ΦΕ111/ἑ5 
Cl=(PlePUe(1.0P31)/12.0eP30ecM) 
۲٤٢١ 
| С2з(Р1#Р2еЕН}/ЕМ 
[δι С3э{Р1әРбеЕМ]/(2.берб) 
CheCle(Clle(C2+2.0eC 31) 
C5-Cle(C11e(C2563)] 


40 
41 
42 
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| STORE Τχ(ὶ 1) FROM 2nd RUN 


— - 

















TCC = EFFECTIVENESS 





CHECK CONVERGENCE 
PAUSE 11 = SOLUTION CONVERSES 
PAUSE 12 = SOLUTION DIVERGES 


REVERSAL CONDITION 


CHECK DIFFERENCE IN EFFECTIVENESS IN 
THE LAST TWC PASSES 


STORE Тх(1,1) + EFFECTIV=NESS AND 


ERROR FOR THE lat RUN 





INCREASE NO. OF SUBDIVISIONS AND REPEAT 
let RUN 


EXTRAPOLATE EFFECTIVENESS AND STORE 


EE —— - 


SUBSCRIPTS FOR ERROR AND EFFECTIVENESS 


URC(I,J) = Te(4,4) + UCC(I,L) = Tare g) 


FOR CONDUCTION CASE 


CONSTANTS FOR CONDUCTICN CASE 








w 5 --- ma سسس ا‎ v-—  —  ———TMH—. Ó—————— аа -y 


53 
ο 


60 
6) 


62 


64 
66 
67 
68 
69 
70 


71 


12 


13 


Tu 


15 





88 
881 


882 


——_ GOTO _ 


3 


--ε м. ce |. BALA — 


01212.0°C1)/C4 

022(CI*(CI11o(2.DeC3-C2)) )/Ch 

03=1C3*C11)/Cu 

Du=(1.0-C1)/C)1 | 
ао CONSTANTS FOR CONDUCTION CASE 

D6=12.0+C11/05 

OT=(Cl+(Clls(>(C3-C2))1/C5 
08=1C3°C111/C5 P | | 


00 67 J2),NH 
L=) | | 
UHC(|,J+))=zDóeTH(1,J)-OTeUHC(),J)+#DB 4 (UFC(2,J)*+UHC(2,J+1)) | 
{ 
i 


M= | ! 
HOT SIDE WITH CONDUCTION | 


TH(2, J) Ob eTHUI,J)*05e(UEFCCUIT; JJ*€UHCC T, Je 2)) 

DOOL 1=2,NMM 
ОНнС!(),)+1)=р}еТН{),4))-Ог»ОНС{ 1 ,Ш))+05е(0НС{1—1!,))+0НС{)-},0)+1}+ 
онС{1+),.))+ОНС{ жт, Је Т 

THU) 1, J)ZOUe THUE,J32*D5e(CUHCCE, J)«UHCCT, Je 1 2) 

UHCINM, J* 1)06 * THULNM, JI) -CZ eLHCIUNE,.N1 *DBe (UHCUNMPM, JJ *UHCENMM, Je 131) ; 


ITERATION ON SIDE OF Cmax 
TWO PASSES PER COLUMN 


THUNMP, J3 QU e THUNM, JI «Doe (UHC(NM, JIJ*€ UHC(NM,.J€ 11) 

|= +1 

YFC3-L) 67,67,61 

CONT INUE 

DO 69 Ks),NM 

UCCIK;, là -UHC( NMP-K , NHP ) 

COLO SIDE WI TH CONOUCT ION 

DO 75 L=1,NC 

Ј= 1 

UCCUI Le) =F OeTCieLI-F7eUCC(),LI+FHe (UCCI2,LIFUCCI2,Lt+I1)1 | 
TCÍ 2; UL ) FU e TCCL T L) *F5e(UCCCIT;U)*UCCEDP Le T 1) 

DO 73 Kz2,NMM 

ОССК, Е ж+ 1) = Е Ie TCUK ,LU-F2*UCCI(K,U) 4F 5e CUCCEK- T,L 1) *UCCCK- 1,0810) 
JUCCUK * IU) *UCCUK* 1,0 1 3) | 
TCUK* ) ,L)SFUeTCUK, LO) *FSe CUCCUK LU *«UCCUK,Le))) | 
UCCUNM,.L* 1)sFÓOeTCUNM, L)-FZeUCC(NM,.L) *F Be LUCCUNMM,.L) *UCCUNMM,L*1)) « 
YC(NMP,.U)SZFUSTCUNM,U) *FS5e (LUCCUNM ,L)FUCCIUNM, Le 1)) 

у= +1 

ırt3-4) 75,75, 71 

CONTINUE 

G1=0.0 

DO 76 Ј=1, мн 

Gl=TH(NMP,¿J1+6G1 

ТНО=СТИЕН 

G220.0 ` 

00 7/8 L=1,NC 


G2=TCINMP,L)+62 
TCO(M1-G2/EC | u 


ITERATION ON SIDE OF Cain 
TWO PASSES PER COLUMN 


ERR(M)=).0-(P)eTCO(M))/( 1.0-TH0) 
)FIERRIM)) 80,80,81 
ERRIM)=-ERRIM) 

IFUERR(M)-4.0E-6) 90,90,E2 
[F(90-M) 83,83,85 





IFCERRIM- 1)-ERREMI) вц,54,87/ 


PAUSE 22 
IFISENSE SWITCH 1} 842,53 


| 
| 
۱۴٢۱١۶۸۶١ SWITCH 2) 88, 884 | 
M=M+ 1 
00 86 1=1,NM 
GO 20 60 
4 


ОНСС, 1) ZUCCENMP- T,NCP ) 
PAUSE 21 


GO TO 1۱1 

ТЕ( ІС-)) 881,881,882 
EF(IA,181=0.D 
ЕК(ТА, 18) =0.0 
50 TO 97 : EF 
FORMATE IIKSOL DIV FOR13, 18X,6E9.2,F8.2,21F)0.6,69.2)) 

WRITE OUTPUT TAPE 2,882,1C, 45, PO ,FF{ TA, IB) ,FROIA,I BI» TCOUM) pe RREM) 





ο 





| 


885 
89 
891 


90 
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IFLSENSE SWITCH 4) 


В 


106 IC-2 


FORMATIUFB.2,E8.2, 68.72, 2X, 135,6 10.6,69.2] 


PRINT 884 ,P1,P2,P3, Pu ,P5,P6, IC, TCOLM1,FRRIM] 


PAUSE 23 

IFUSENSE SWITCH 1) 885,53 
IF(SENSE SWITCH 2) 114,8€ 
۱۴۶۲1٢۰٠ 891.891,93 
٤٤٢) ))۱ 

6 8 ۱ 

66.10 182 4205 
TI=TCOUM-1)-TCOUM) 

IFUT H1 91,92,92 

Ti=-11 

IFUTI-4.0E-61 89,89,82 


EFF=TCO[MI=*(1.0+REI-EF(1A, IF 1+RE 
82=(1R1-EFF)+100.0)/R 1 


1,82 


IFISENSE SWITCH 31 95,102 
FORMAT(F8.2+E9.2+£8.41 
PRINT 95,P4,P5,82 
60 TO 102 . 

00 97% 1=),мм 

00 972 J=1,NHP 

UHC (1, J)=UH( 10,3) 

DO 97% (-1,МСР 
UCC(l1,L1=UC(I,L) 
131 132 
FORMATU/4EB.2+E8.2+F8.213X+13/(115F7.4)) 


WRITE OUTPUT TAPE 32131,P1,P2,P4,PU,P5,P6,I0, 
HUHCCJ, VI UCCUNMP- Ju NCP) , Je T, NH) 


Ρό-Ρό-ΡόΙΝ 

1865 [89 ἱ 

IF(P6-P6FI) 103,103,51 
P6=S6 

PS=PSePSFAtPSIN 
lA=1A+) 

IFtPSFI-PSI 105,50,50 
Р5-55 

IF(1C-1) 106, 106, 108 


NM=NM2 

NH=NH2 

NC=NC2 

DO 107 K=1,NM 
UH(K, 1) 2UH24K ) 

GO TO 9 

IF(PUFI-=P41 113, 113,109 
PuzPue«Pu IN 


NMzNM-NM| 
NH=NH-NHI 


NC-NC-NCI 
00 111 K=1,NM 
UHIK, I1=UHlIK) 


Pr‏ امت سے وجب ۲٢۰۹۰‏ و 


PAUSE 50 

IFISENSE SWITCH 11 llel 
DO 115 L=1,NM 

UHLL, I1=UH2LL 1 

σο το 5 





ΕΝΟ 
ENO 


FORMAT(32X,EB.2,FB.2,2(F 10.6, E9.2 1, F 10.6,F B. 31 
WRITE OUTPUT TAPE 2,95 ,P5,P65, EFCITA,IRT,ERCIA,S IBI S TCOULUV)  ERRIMP IT, EFF 
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